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Many non-equilibrium systems display dynamic phase transitions from active to absorbing 
states, where fluctuations cease entirely. Based on a field theory representation of the master 
equation, the critical behavior can be analyzed by means of the renormalization group. The 
resulting universality classes for single-species systems are reviewed here. Generically, the 
critical exponents are those of directed percolation (Reggeon field theory), with critical di- 
mension d c = 4. Yet local particle number parity conservation in even-offspring branching 
and annihilating random walks implies an inactive phase (emerging below d' c ~ 4/3) that is 
characterized by the power laws of the pair annihilation reaction, and leads to different critical 
exponents at the transition. For local processes without memory, the pair contact process with 
diffusion represents the only other non-trivial universality class. The consistent treatment of 
restricted site occupations and quenched random reaction rates are important open issues. 

PACS: 64.60.Ak, 05.40.-a, 82.20.-w 

1 Introduction: Active to absorbing state phase transitions 

Among the prevalent goals in current statistical mechanics is the understanding and character- 
ization of non-equilibrium steady states. This program is hindered by the general absence of 
an effective free-energy function that would allow a straightforward classification in terms of 
symmetries and interactions. One might hope, however, that the task becomes more feasible 
near continuous phase transitions separating different non-equilibrium steady states. For in anal- 
ogy with equilibrium critical points, one would expect the properties near non-equilibrium phase 
transitions as well to be universal, i.e., independent of the detailed microscopic dynamical rules 
and the initial conditions. Rather, the emerging power laws and scaling functions describing the 
long-wavelength, long-time limit should hopefully be characterized by not too many distinct uni- 
versality classes. Obviously, provided we can cast the problem at hand into a form amenable to 
field-theoretic methods, the renormalization group (RG) provides a very powerful tool for such 
investigations. Indeed, during the past twenty-five years or so it has been successfully applied to 
a variety of non-equilibrium processes. Among the lessons we have learned is that critical phe- 
nomena or generic scale invariance far from thermal equilibrium, where the detailed-balance con- 
straints do not apply, are considerably richer than equilibrium statics, or even near-equilibrium 
dynamics. In fact, intuitions drawn from the latter may often be quite deceptive. 
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A special class of genuine non-equilibrium phase transitions separate 'active' from 'inactive, 
absorbing' stationary states where any stochastic fluctuations cease entirely [1,2]. These occur 
in a large variety of systems, e.g., in chemical reactions involving an inert state that does not 
release the reactants A anymore. We may also consider stochastic population dynamics, com- 
bining, say, diffusive migration with asexual reproduction A — > 2A (with rate a), spontaneous 
death A —> (rate fi), and lethal competition 2A — > A (rate A). In the inactive state, where 
no population members A are left, all processes terminate. Similar effective dynamics may be 
used to model non-equilibrium physical systems, such as the domain-wall kinetics in Ising chains 
with competing Glauber (spin flip) and Kawasaki (spin exchange) dynamics [3]. Here, spin flips 
TTII^TTTi an d TTIT ~~ > T' T mav be viewed as domain wall (A) hopping and pair annihila- 
tion 2 A — > 0, respectively, whereas a spin exchange TT-i-i — * TITI represents a branching process 
A — > 3 A in domain wall language. Notice that the paramagnetic and ferromagnetically ordered 
phases map onto the active and inactive 'particle' states, the latter rendered absorbing if the spin 
flip rates are computed at zero temperature, allowing no energy increase. 

The simplest mathematical description for such processes uses a kinetic rate equation in terms 
of the time-dependent average 'particle' density n(t), which for the above systems reads 

d t n{t) = (a-n)n(t) - Xn{tf . (1) 

Obviously, this yields both an inactive and an active phase, as for a < \i we have n(t — > oo) — ► 0, 
whereas for a > fi the particle density saturates at n s = (a — fi)/ A. The explicit solution n(t) = 
n n s / [no + (n s — no)e^ _<7 )*] shows that either stationary state is approached exponentially in 
time. The two phases are separated by a continuous dynamic transition at the critical point u = fi, 
where the temporal decay becomes algebraic, n(t) = uq/{1 + no At). However, eq. (1) entails 
a mean-field type of approximation, as we have neglected particle correlations on the right- 
hand side, and effectively factored a two-point correlation function, namely the joint probability 
of finding two particles at the same position. A more detailed treatment therefore requires a 
systematic incorporation of spatio-temporal fluctuations, and the ensuing particle correlations. 



2 Langevin description, Reggeon field theory, and directed percolation (DP) 

We may try a phenomenological incorporation of fluctuation effects using a Langevin-type ap- 
proach. To this end, we assume diffusive particle transport, and model our 'chemical' system 
through a non-linear stochastic differential equation with reaction functional r[n], 

d t n(x, t) = DV 2 n(x, t) - r[n] (x, t) + C(x, t) . (2) 

The stochastic variable £(x, t) incorporates the reaction noise, which we take to have zero mean, 
and to be local in space-time ('white' in Fourier space), with a density-dependent correlator c[n]: 

(0=0, (C(x,i)C(x',0>=c[n]5(x-x')^-iO- (3) 

Here, the last expression means that whenever a trajectory average is taken, c[n] needs to be 
factored in. In the spirit of Landau theory, we may expand the functionals r[n] and c[n] near the 
inactive phase (n <c 1). In the absence of spontaneous particle production, both must vanish at 
n = 0, which is the condition for an absorbing state. Consequently, an active to absorbing state 
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phase transition should be generically described by eqs. (2), (3) with r[n] « Dm + \n 2 and 
c[n] m 2an [4]. Upon neglecting fluctuations, one then recovers eq. (1) with r = (p, — a)/D. 

By means of standard techniques, a stochastic differential equation of the form (2) with noise 
correlation (3) can be represented as a functional integral [5]. Essentially, one starts with the 
Gaussian noise probability distribution W[Q oc exp [— J((^ 2 /2c)d d x dt] , multiplies with 1 = 
/ D M Y[^. t 5{dtn-D\/ 2 n+r-C) = J D[in]D[n]exp [- J h(d t n - DV 2 n + r - ()d d x dt], 
thus introducing the auxiliary fields h, and then integrates out the stochastic noise (. Upon 
applying a forward discretization, the ensuing functional determinant vanishes, and one ar- 
rives at the probability distribution P[n] oc J D[ih]exp(—S[n,n]), with the response func- 
tional S[n, n] — J h (d t n — DV 2 n + r[n] - c[n]n/2) d d xdt. After rescaling n = (a/X) 1 / 2 ^, 
n = (A/ct) 1 / 2 ^, the expanded Langevin equation (2), (3) maps to 'Reggeon' field theory [6] 

d d xdt, (4) 

with u = (erA/2) 1 / 2 . Notice the invariance of this action with respect to 'rapidity inversion' 

0(x, t) -> -0(x, -t), 4>(x, t) -> -0(x, -t). 

The field theory (4) should capture the generic critical behavior for non-equilibrium phase 
transitions between active and absorbing states, occurring at r = 0. Quite remarkably, the very 
same action is obtained for the threshold pair correlation function [7] in the purely geometric 
problem of directed percolation (DP) [8]. Power counting reveals d c = 4 as the upper critical 
dimension; hence, the critical exponents as predicted by mean-field theory acquire logarithmic 
corrections at d c , and are shifted to different values by the infrared-singular fluctuations in d < 4 
dimensions. By means of the standard perturbational loop expansion in terms of the diffusion 
propagator and the vertices oc u, and the application of the RG, the critical exponents can be 
computed systematically and in a controlled manner in a dimensional expansion with respect to 
e = 4 — d. The one-loop results, to first order in e, as well as reliable values from Monte Carlo 
simulations in one and two dimensions [2] are listed in Tab. 1 . Moreover, as a consequence of 
rapidity invariance there are only three independent scaling exponents, namely the anomalous 
field dimension r/, the correlation length exponent v, and the dynamic critical exponent z. All 
other exponents are then fixed by scaling relations, such as j3 = v(d + z — 2 + ry)/2 = zva for 
the order parameter exponent (3 and the critical density decay exponent a, respectively. 



DP exponents 


d= 1 


d = 2 


d = 4 - e, O(e) 


n s ~ \r\ p 
£~ \r\~" 
£ c ~e~ \r\-* v 
n c (t) - t~ a 


(3 w 0.2765 
v w 1.100 
z « 1.576 
a w 0.160 


13 w 0.584 
v « 0.735 
z w 1.73 
a w 0.46 


(3 = 1 - e/6 
= l/2 + e/16 
z = 2 - e/12 
a = 1 - e/4 



Tab. 1 . Critical exponents for the saturation density (order parameter) n s , correlation length £, characteristic 
time scale £ c , and critical density decay n c (t) for the universality class of directed percolation (DP). 
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3 From the master equation to stochastic field theory 



The above phenomenological approach constitutes a natural extension of the mean-field rate 
equations to a stochastic partial differential equation. However, it does presume (i) that such a 
Langevin-type representation is in fact possible, and (ii) it is fundamentally based on conjectures 
on the noise correlator. Away from thermal equilibrium, there is no analog of the Einstein relation 
to constrain the form and structure of the noise correlations. Indeed, it has emerged that the latter 
may quite profoundly affect the scaling behavior of non-equilibrium systems, rendering point 
(ii) an entirely non-trivial issue. Moreover, as we shall see, even assertion (i) turns out to be 
relevant in certain important model systems. It is therefore of fundamental importance to be 
able to construct a long-wavelength or field theory representation of stochastic processes that 
starts directly from their microscopic definition in terms of a classical master equation, without 
recourse to any serious additional assumptions or approximations. 

Fortunately, for reaction-diffusion systems there exists indeed a standard route from the mas- 
ter equation to an effective 'Hamiltonian' (more precisely, the Liouville time evolution operator), 
and therefrom immediately to a field theory action [9]. The key point is that all possible con- 
figurations here can be labeled by specifying the occupation numbers rii of, say, the sites of a 
rf-dimensional lattice; we shall henceforth assume that there are no site occupation restrictions, 
i.e., rii = 0, 1, 2, . . . The master equation then addresses the time evolution of the configurational 
probability P({n,}; t). For example, the corresponding contribution from the binary coagulation 
process 2A — ► A at site i reads d t P{rii; t)\ \ — A [(m + l)rijP(nj + l;i) — nj(rij — l)P(n,;i)]. 
This sole dependence on the integer variables {rii} calls for a second-quantized bosonic oper- 
ator representation with the standard commutation relations [a*, at] = Sij and the empty state 

|0) such that a^O) = 0. We then define the Fock states via |{nj) = n»( a i) n< l°) (notice 
that the normalization is different from standard many-particle quantum mechanics), and thence 
construct the formal state vector |$(t)) = X){n } P{{ n i]\ t) \{ n i})- The master equation now 
imposes a linear time evolution that can be written as an imaginary-time 'Schrodinger' equation 
9 t | <!>(£)) = —H |$(£)), with a generally non-Hermitian stochastic 'Hamiltonian' H ({at}, {a^}); 
e.g., for the on-site coagulation reaction one finds explicitly H\ j = — A(l — a\)a\af. 

Our goal is to evaluate time-dependent statistical averages for observables F, necessarily 
mere functions of the occupation numbers as well, whence (F(t)) — J2{ ni } F({ n i})P({ n i}'i 
Straightforward algebra utilizing the identity [e a , a*] — e° shows that this average can be cast 
into a 'matrix element' (F(t)) = (V\F({ai})\$(t)) = {P\F({a t }) e - Ht |$(0)} with the state 
vector |3>(i)) and the projector state (V\ = (0 J\ i e ai , with (V\0) = 1. For example, proba- 
bility conservation implies 1 = (V\e~ Ht |3>(0)), i.e., for infinitesimal times (V\$>(0)) = 1 and 
(V\H = 0, which is satisfied if iJ({l},{aj}) = 0. Notice furthermore that commuting the 
factor * 0i through all the other operators has the effect of shifting a\ — > 1 + a\ everywhere. 

As a final step, we employ the coherent-state path integrals familiar from quantum many- 
particle systems [10], and perform the continuum limit to arrive at the desired field theory. For 
our earlier population dynamics example with diffusive motion and the reactions A — > 2A, 
A — > 0, and 2 A — > A, this results in the action (omitting contributions from the initial state) 



S$, ip] = J j>(dt- DV 2 )ip + a(l - tp)W - n(l - ip)ip - A(l - i>)^ 2 



d d xdt. (5) 
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It is important to emphasize again that the single approximation used here was the continuum 
limit; specifically, no assumptions whatsoever on the stochastic noise were invoked. The mean- 
field rate equation (1) is recovered by inserting the solution r/> = 1 to the stationarity condition 
SS/St/j = into SS/8tp = 0, with n(t) — 2(ip(x,t)). However, higher moments of the field 
xp cannot be directly identified with the corresponding density correlations. Very remarkably, 
though, after shifting ip(x, t) = 1 + ip(x, t), such that {ip) — 0, and appropriate field rescaling 
one arrives again at the field theory (4) with r = (ji — cr)/D and u = {2a\) 1 / 2 , with the 
additional vertex X(f> 2 (j) 2 . Yet, the coupling A has scaling dimension n 2 ~ d , where n denotes a 
momentum scale, whereas u becomes marginal at d c = 4. Thus, at least in the vicinity of the 
upper critical dimension, A constitutes an irrelevant variable in the RG sense, and in general the 
ratio \/u ~ KT d l 2 is expected to scale to zero asymptotically. Consequently, for processes in the 
directed percolation universality class the microscopic master equation and the coarse-grained 
Langevin description both lead to the effective action of Reggeon field theory. 



4 Diffusion-limited annihilation 

Let us now investigate the fc-th order annihilation reaction kA — ► [11]. The corresponding 
mean-field rate equation reads d t n(t) = —\n(t) k . For simple radioactive decay (k = 1), it is 
solved by the familiar exponential n(t) = n n e~ Al *, whereas one obtains power laws for k > 2, 

namely n(i) = \n]^ k + (k — 1) A t\ ^ \ln order to consistently include fluctuations in the 
latter case, we start out from the master equation once again, which leads to the action 

S$ t $\ = J i> (d t - DV 2 ) ip - \(l - ip k )ip k d d xdt. (6) 

After shifting tp(x,t) = 1 + ip(x, t), it becomes clear that this field theory actually has no simple 
Langevin representation. For in order to interpret tp as the corresponding noise auxiliary field, it 
should appear quadratically in the action only, and with negative prefactor. Thus, the Langevin 
equation derived from the action (6) for the pair annihilation process would imply unphysical 
'imaginary' noise with c[n] = —2\n 2 . 

Analyzing the field theory (6) further, we see that the diffusion propagator does not become 
renormalized at all, implying that rj = and z = 2 to all orders in the perturbation expansion. 
The critical dimension of the annihilation vertex is found to be d c (k) = 2/(k — 1), allowing 
for the possibility of non-trivial scaling behavior in low physical dimensions only for the pair 
and triplet processes. The simple structure of the action permits summing the entire perturbation 
series for the vertex renormalization by means of a Bethe-Salpeter equation that in Fourier space 
reduces to a geometric series [11]. For pair annihilation (k = 2), this yields the following 
asymptotic behavior for the particle density: n(t) oc t^ 1 for d > 2, n(t) oc t^ 1 Int at d c = 2, 
and n(t) oc t~ d l 2 for d < 2. The slower decay for d < 2 originates in the fast mutual annihilation 
of any close-by reactants; after some time has elapsed, this leaves only well-separated particles. 
The annihilation dynamics thus produces a«f/-correlations, mimicking effective repulsion (which 
is also the physical interpretation of the negative sign in c[n]). Similarly, for triplet annihilation 
(k = 3) the density decays as n(t) oc i -1 / 2 for d > 1, with mere logarithmic corrections 

n(t) oc (t- 1 hit) 1/2 at 4 = 1. 
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5 Branching and annihilating random walks (BARW) 

In order to allow again for a genuine phase transition, we combine the annihilation kA — > 
(k > 2) with branching processes A — > (m + 1) A. The associated rate equation reads dt n(t) = 

-Xn{t) k + an(t), with the solution n(t) = n s / (l + [(n s /n ) k - 1 - l] e -( fe - 1 ) <rt ) 1/(fe_1) . 
Mean-field theory thus predicts the density to approach the saturation value n s = (a/X) 1 ^ k ~ 1 ^ 
as t — > oo for any positive branching rate. Above the critical dimension rf c (fc) = 2/(k — 1) 
therefore, the system only has an active phase; a c = represents a degenerate 'critical' point, 
with scaling exponents essentially determined by the pure annihilation model: a = (3 = l/(k — 
1), v = 1/2, and z = 2. However, Monte Carlo simulations for these branching and annihilating 
random walks (BARW) revealed a much richer picture, in low dimensions clearly distinguishing 
between the cases of odd and even number of offspring m [3, 12]: For k — 2, d < 2, and m odd, 
a transition to an inactive, absorbing phase is found, characterized by the DP critical exponents. 
On the other hand, for even offspring number a phase transition in a novel universality class with 
a s=s 0.27, (3 w 0.92, v « 1.6, and z w 1.75 emerges in one dimension. 

The aforementioned mapping to a stochastic field theory, combined with RG methods, suc- 
ceeded to elucidate the physics behind those remarkable findings [13]. The action for the most 
interesting pair annihilation case becomes 



S[i>, *!>]=[ $ (dt - D V 2 ) ip - A (1 - ip 2 ) tp 2 + a (1 - m ) ^ d d x dt . 



(7) 



Upon combining the reactions A — > (m + 1)A and 2A — > 0, one notices immediately that the 
loop diagrams generate the lower-order branching processes A — > (m — 1)A, A — > (m — 3) A . . . 
Moreover, the one-loop RG eigenvalue y a = 2 — m(m + l)/2 (computed at the annihilation 
fixed point) shows that the reactions with smallest m are the most relevant. For odd m, we see 
the generic situation is given by m — 1, i.e., A — ► 2A, supplemented with the spontaneous decay 
^4^0. After a first coarse-graining step, this latter process (with rate /i) must be included in 
the effective model, which hence becomes identical with action (5). We are thus led to Reggeon 
field theory (4) describing the DP universality class, provided the induced decay processes are 
sufficiently strong to render a c > 0. Yet for d > 2 the renormalized mass term <tr — /i^ 
remains positive, which leaves us with merely the active phase captured by mean-field theory. For 
d < 2, however, the involved fluctuation integrals are infrared-divergent, thus indeed allowing 
the induced decays to overcome the branching processes to produce a phase transition. As the 
DP upper critical dimension is d c = 4, the scaling exponents display an unusual discontinuity 
at d = 2, jumping from the non-trivial two-dimensional DP to the mean-field values (for any 
d > 2) as a consequence of the vanishing critical branching rate [13]. 

It is now obvious why the case of even offspring number m is fundamentally different: Here, 
the most relevant branching process is A — > 3A, and spontaneous particle death with associ- 
ated exponential decay is not generated under the coarse-grained dynamics, which precludes the 
above mechanism for producing an inactive phase with exponential particle decay. This impor- 
tant distinction from the odd-m case can be traced to a microscopic local conservation law, for 
the reactions 2 A — > and A — ► 3A, A — ► 5 A . . . always destroy or produce an even number 
of reactants, and consequently preserve the particle number parity. Formally, this is reflected in 
the invariance of the action (7) under the combined inversions ip — ► — tp, -0 — > — ip. As we saw 
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earlier, the branching rate o certainly constitutes a relevant variable near the critical dimension 
d c = 2. Therefore the phase transition can only occur at o c = 0, and for any o > there exists 
only an active phase, described by mean-field theory. In two dimensions one readily predicts the 
following logarithmic corrections: £(<r) oc er -1 / 2 ln(l/cr), and n(o) oc o [ln(l/cr)]~ 2 . How- 
ever, setting m = 2 in the one-loop value for the RG eigenvalue y a , we notice that the branching 
vertex becomes irrelevant for d < 4/3. More information on the low-dimensional behavior can 
be gained through a one-loop analysis at fixed dimension, albeit uncontrolled [13]. The ensu- 
ing RG flow equations for the renormalized, dimensionless branching rate or = <j/Dk 2 , and 
annihilation rate \ R = C d \/DK 2 ~ d , with C d = Y{2 - d/2)/2 d - 1 n d / 2 read (form = 2): 



da R 



2 3A * 



(1 + o R ) 2 ~ d / 2 



2-d- 



(I + Or) 2 ^/ 2 



(8) 



The effective coupling is then identified as g = Xr/(1 + or) 2 ^ 4 / 2 , which approaches the an- 
nihilation fixed point g* = 2 — d as or — > 0, while for or — > oo the flow tends towards the 
Gaussian fixed point g* = describing the active state. The separatrix between the two phases is 
given by the unstable RG fixed point g* = 4/(10 — 3d), which enters the physical regime below 
the novel critical dimension d' c w 4/3. For d < d' c , this describes a dynamic phase transition 
with o c > 0; the emerging inactive phase is characterized by a vanishing branching rate, and 
thus by the algebraic scaling behavior of the pure pair annihilation model, n(t) oc t~ d / 2 . This 
very different character of the inactive phase also explains why the transition cannot possibly fall 
into the DP category. The aforementioned fixed-dimension RG analysis yields the rather crude 
values v « 3/(10 — 3d), z w 2, and (3 w 4/(10 — 3d) for this parity-conserving (PC) universality 
class. Moreover, the absence of any mean-field ounterpart for this transition precludes a sound 
derivation of 'hyperscaling' relations such as [3 — zva. 

Similar arguments suggest for triplet annihilation 3 A — > combined with branching A — > 
(m+l)A that DP behavior with o c > should ensue for m mod 3 = 1, 2, as then the processes 
A — > 0, A — > 2A, and 2A — > A are dynamically generated. For m = 3, 6, . . ., on the other hand, 
because of d c (k = 3) = 1 one expects logarithmic corrections only in one-dimensional systems. 



6 Concluding remarks and outlook 



In this brief overview, I have outlined how non-linear stochastic processes can be represented by 
field theory actions, allowing for a thorough analysis and classification by means of the renor- 
malization group. Systems with a single 'particle' species displaying a non-equilibrium phase 
transition from an active to an inactive, absorbing state are generically captured by the DP univer- 
sality class. The second prominent example, applicable when additional symmetries (degenerate 
absorbing states) are present, is that of branching and annihilating random walks with even off- 
spring number (PC). In fact, there appears to be only one additional possibility for non-trivial 
scaling behavior in 'bosonic' systems with multiple site occupation, namely the combination of 
binary reactions 2A — > and 2A — > (m + 2)A (annihilation and fission; pair contact process 
with diffusion, PCPD) [14]. These reactions subsequently generate 2 A — ► mA, (m — 2) A, . . ., 
2(m + 1)A, . . ., thus producing infinitely many couplings with identical scaling dimensions, 
which renders the field theory non-renormalizable. Yet, the master equation immediately leads 
to the exact evolution equation d t n(t) = (o — X)(m(m — 1)), which establishes the existence 
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of a phase transition when the particle production and annihilation rates precisely balance. In the 
inactive, absorbing phase, obviously the power laws of the pure annihilation (m even) or coagu- 
lation (to odd) models are recovered. In the active phase, the particle density diverges in a finite 
time, and the transition itself is discontinuous. These irregular features are 'cured' by restricting 
the site occupation numbers, whereupon the phase transition becomes continuous, with critical 
exponents a « 0.13, j3 « 0.58, z w 1.83 that belong to neither the DP nor PC universality 
classes. Yet different exponents were recently found for the triplet reactions 3A — > 2A and 
3 A — > A A [16]. It is however unclear at present how to generally and systematically incorporate 
such site restrictions into a field theory description in a tractable manner [17]. Another open 
problem concerns the obviously very relevant issue of quenched disorder in the reaction rates. 
For example, the RG investigation for DP with random percolation threshold yields run-away 
flows to diverging couplings, with an as yet unclear interpretation [18]. 

Lastly, the rich variety of multi-component systems is only partially understood. Only DP 
processes with arbitrarily many particle species have been fully classified, namely any non-linear 
coupling leads directly to the DP universality class again [19], whereas novel multi-critical be- 
havior may ensue for hierarchical systems comprising spontaneous transformations A — > B, 
B — > C . . . [20, 19]. Yet already in the simple binary process A + B — > 0, the initial condi- 
tions play a vital role. For equal A and B densities, the fact that ua — tib remains conserved 
under the dynamics leads to segregation and reaction-front dominated kinetics, with shifted up- 
per critical dimension d c = 4 and asymptotic decay i%a oc jib oc t~ d / 4 for lower dimensions 
[21]. On the other hand, whenever there is a majority species, the minority decays to zero expo- 
nentially, although this may be masked by very long crossover times. The (bosonic) TV-species 
BARW generalization Ai + Ai — > 0, A, — > Ai + 2Aj turns out to be exactly analyzable, with 
a degenerate phase transition only at vanishing branching rate, and critical exponents a = d/2, 
(3 = 1, v = 1/d, z = 2 [13]. The situation for N = 1 is thus qualitatively different from 
all multi-component cases. Moreover, in certain cases the difference in diffusion constants is 
a relevant control parameter [22]. When there are 'passive agents, the corresponding variables 
are readily integrated out, leaving however interactions of the remaining degrees of freedom that 
are non-local in time. A classic example for such memory effects is dynamic percolation [23], 
which at criticality yields stationary isotropic percolation clusters. A recent argument claims 
this to govern the generic universality class for non-equilibrium phase transitions with infinitely 
many absorbing states [24]. Many intriguing problems are still open, and RG techniques will be 
invaluable tools for the further analysis of scale-invariant non-equilibrium systems. 
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